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Abstract 

Axially symmetric monopole anti-monopole dipole solutions to the second order equations of a simple 
SU(2) Yang-Mills-Higgs model featuring a quartic Skyrme-like term are constructed numerically. The effect 
of varying the Skyrme coupling constant on these solutions is studied in some detail. 
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1 Introduction 



The SU(2) Georgi-Glashow model in the Bogomornyi-Prasad-Sommcrficld (BPS) limit supports monopoles [1, 2] 
which are solutions of the first order self-duality equations [3, 4]. Away from the BPS limit, when new gauge 
invariant and positive definite terms are added, the resulting monopoles are described by the solutions to the 
second order Euler-Lagrange equations, and not to the first order self-duality equations. Once these terms are 
introduced to the model, the BPS topological bound cannot be saturated. 

BPS and non-BPS monopoles differ in two remarkable respects. First, the BPS multimonopoles can be 
constructed analytically [6, 7, 8, 9] while the non-BPS monopoles, e.g. when the Higgs potential is present [1, 2], 
can only be constructed numerically. Secondly, and perhaps physically more interestingly, BPS monopoles do not 
interact while non-BPS monopoles interact. In the presence of a Higgs potential this interaction is known to be 
repulsive [10, 11] and has been verified to be so numerically [12], while in the presence of Skyrme like terms, higher 
order in both the Yang-Mills (YM) curvature and the Higgs covariant derivatives, this interaction can be both 
repulsive and attractive [13]. In a particularly simple such (Skyrme like) model, this interaction was found [14] to 
be strictly attractive, and moreover it was found [14], rather unexpectedly, that the lowest energy bound states 
were the axially symmetric ones and not those with Platonic symmetries. (It was unexpected since this feature 
contrasts with that for Skyrmion bound states [15].) 

All the above monopole solutions discussed are stable relative to the topological lower bound whether they 
saturate this bound, as for the BPS monopoles, or not, as for non-BPS ones. There is however another class of 
non-selfdual solutions to the second order Euler-Lagrange equations which are not stable and represent states 
of monopoles and anti-monopoles in equilibrium. The existence of such solutions was first proved by Taubes [17] 
for the model featuring no Higgs potential (and of course no higher order terms in the curvature and covariant 
derivative), namely for the model which supports BPS multimonopoles. Such a non-BPS solution, namely an 
unstable solution of the second order equations, was first constructed for this system with SU(3) gauge group and 
subject to spherical symmetry by Burzlaff [18]. More recently Ioannidou and Sutcliffe [19] emoplyed a harmonic 
map Anstaz to construct such spherically symmetric solutions to the same (BPS) system with gauge groups 
SU(3), SU(4) and SU(N). Using results on sigma model instantons, these authors [19] also argued that the zero 
charge solutions they constructed described monopole anti-imonopole pairs. 

A direct approach to constructing zero charge monopole anti-monopole pairs for the SU (2) BPS model was 
used sometime ago by Ruber [20]. This was the numerical construction of axially symmetric solutions with 
suitable boundary value conditions. More recently Klcihaus and Kunz [21] constructed this zero charge solution 
for the full Georgi-Glashow model featuring a Higgs potential, and they studied the effect of the Higgs potential 
in detail. To date, no such study has been reported in the literature pertaining to the model featuring higher 
order Skyrme like terms. In the background of the above described scenario it is pertinent to carry out such a 
study. 

This is the aim of the present work. We will consider the zero charge axially symmetric monopole anti- 
monopole solutions as in [20, 21], for the simple skyrmed Higgs model studied in [14] whose axially symmetric 
charge-2 monopoles are mutually attractive. This contrasts with the monopole anti-monopole solutions studied 
in [21] for the model whose charge-2 monopoles are mutually repulsive, which makes the comparison of our results 
with those of [21] interesting. In addition to constructing the vorticity-1 monopole anti-monopole solutions, as in 
[20] and [21], but now for the skyrmed model here, we also construct the corresponding vorticity-2 solutions. 



2 Skyrmed SU(2) Yang-Mills-Higgs Model 



The static energy of the simplified Skyrme like model considered is 




d?r (1) 



with field strength tensor of the su(2) gauge potential = \r a A a ^ 

Fp V = d^Av - d v A^ + ig[A^,A v \ , 
and covariant derivative of the Higgs field <I> = r a (j) a in the adjoint representation 

= d^ + ig[A^] , 



(3) 



(2) 
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and g denotes the gauge coupling constant, k the coupling strength of the quartic Skyrme like Higgs kinetic term, 
A the strength of the Higgs potential and -q the vacuum expectation value of the Higgs field. 
The topological charge Q is the well known quantity 

Q = -t-/ ik I ^ r {FijDk®} d 3 r j (4) 

corresponding to the magnetic charge m — Q/g, and takes integer values that equal the winding number of the 
Higgs field [22]. The latter is encoded with the boundary conditions which yield the value of this integer. 

To construct axially symmetric solutions that describe systems of monopoles and multimonopoles, specific 
boundary conditions must be imposed the Higgs field at infinity. For usual multimonopoles, the Higgs field at 
infinity is described by the vortex number n winding the azimuthal angle p, n times and the polar angle 9 does 
not wind. Zero magnetic charge monopoles on the other hand, namely those we seek to construct, can be achieved 
by requiring that in the asymptotic Higgs field the polar angle is enhanced by another integer m. This can also be 
achieved automatically by incorporating this integer m in the Ansatz [20, 21] as will be done below. The integral 
(4) can be evaluated for a system with m zeros of the Higgs field (i.e. with m monopole and antimonopole 
centres), and with vorticity n, yielding 

Q = 47TTO7 3 [l-(-in . (5) 

In this paper we will restrict to the charge zero case m = 2 with vorticity n = 1, to carry out our detailed analysis 
of the system, with special attention to the k dependence of the solutions. After that, we will briefly study also 
the case of n = 2 vorticity, again with m = 2. These are both monopole anti-monopole solutions to the second 
order equations carrying Q = 0. 



3 Static axially symmetric Q = Ansatz 

We choose the static, axially symmetric, purely magnetic Ansatz employed in [20] for the monopole-antimonopole 
solution and in [23, 24] for the sphaleron-antisphaleron solution of the Weinberg-Salam model. Here the gauge 
field is parametrized by 

A = , A r = |f , A g = M^), Av = (ff 3T (2,n) + (1 _ H4) T ( 2 .«)) , (6) 

and the Higgs field by 

* = V (*i4 2 ' n) +$ 2 rf n) ) . (7) 
All functions H\, H 2 , H3, i/4,<I?i and $2 depend on (r,8) or equivalently on (p = rs'md,z — rcosd), with the 

An 
Tip 



su(2) matrices Tr 2 ' n \ t^ 2 '™^ and defined in terms of the Pauli matrices ti,t 2 ,t 3 as 



T (?,n) _ gm 2#(cos nip n + sin npT 2 ) + cos 29 T3 , 
T ( 2 ' n ) _ cos 2Q (cos n<p t\ + sin npT 2 ) — sin 29 r 3 , 
r^") = — sin Tip 7~1 -\- COS Tip T2 , (8) 

and for later convenience we define 

cos Tip 7~i -(- sin Tip T 2 . (9) 



r, 



(n) 



P 

Note that the dependence on the vorticity n is encoded through Tr 2 '™-* and Tg 2 ' n \ and of course Tp n \ 
We change to dimcnsionlcss coordinates, Higgs field and coupling parameters by rescaling 

grj g 2 Tj 4 g 2 

respectively. Then this Ansatz leads to the field strength tensor 

Fre = ~{d e H 1 +2rd T H 2 )Tp> , 
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F rip = { (sin 29H X - 2 sin 0i?i(l -H A )-2 sin 6V<9 r iJ 3 ) t^ 2 '^ 

+ (cos 26H! + 2 sin 0i2i# 3 + 2 sin 6rd r H A ) rf ' n) } , 
F ev = {(2 sm29(H 2 -l) + 2 cos 6H 3 - 2 sm9H 2 (l - H 4 ) + 2 sm9d e H 3 )T^ 

+ {2cos29(H 2 - 1) + 2 cos (9(1 - H A ) + 2sin9H 2 H 3 - 2sin0<9 e tf 4 ) rf'" 5 } , (10) 
and the covariant derivative of the Higgs field 

£> r $ = i {(rd r $i + Hi$ 2 ) ^ 2 ' n) + (rd r $2 - r e (2 ' n) } , 

= (9 e $i - 2tf 2 $ 2 ) r^ 2 '") + (a e $ 2 + 2ff 2 $i) 4 2 ' n) , 

= n {(sin 29 -2 sin (9(1 -tf 4 ))$i + (cos 20 + 2 sin (9iJ 3 ) $ 2 } T^ n) . (11) 
The dimensionless energy density then becomes 

+ i T r {0D r $) 2 + 1(^$) 2 + -^(^*) 2 } 

- ^Tr (l[£> r *,£> 9 $]2 + -^W£> r $,£> v $p + - T -L-[ J D e $,£) ¥ ,$] 2 | + A((|<i>| 2 - l) 2 , (12) 
4 (_ r z r 2 sin 6> r 4 sin 9 J 

where $| = ^/<f> 2 + $ 2 denotes the modulus of the Higgs field. 

For a monopole-antimonopole pair we expect a magnetic dipole field for the asymptotic gauge potential. 
The dipole moment C m can be extracted from the gauge field function H 3 , in the gauge where the Higgs field 
approaches asymototically a constant. Like in Ref. [21] we find 



ff 3 = — sin0, (13) 
r 

while all other gauge field functions decay faster. 

4 Numerical Results 

As noted in [21] the Ansatz Eqs. (6), (7) possesses a residual U(l) gauge symmetry. To obtain an unique solution 
we use the gauge fixing condition [21] 

G t = \ (rdrHx - 2dgH 2 ) = . (14) 

The system of partial differential equations is solved numerically subject to the following boundary conditions, 
which respect finite energy and finite energy density conditions as well as regularity and symmetry requirements. 
These boundary conditions are at the origin 

H 1 (0,6)=H 3 (0,6) = , H 2 (0,6) = H 4 (0,e) = l , (15) 

sin20$i(O,(9) + cos2<9$ 2 (O,0) = , <9 r (cos 20$i(O, 9) - sin 2(9$ 2 (0, 9)) = , (16) 

at infinity 

i3i(oo,0) =i? 2 (oo,0) =0 , # 3 (oo,6>) = sin0 , (1 - #4(00,0)) = cos<9 (17) 

$i(oo,0) = 1, $2(^0,60 = 0, (18) 

and on the z-axis 

#i(r,0 = O,7r) = tf 3 (r,0 = O,7r) = d e H 2 (r,6 = 0,n) = d e H A (r, 9 = 0, tt) = , (19) 
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$ 2 (r,0 = 0,tt) = <9o$i(r,6> = 0, tt) = . (20) 

The numerical calculations were performed with the software package CADSOL, based on the Newton- Raphson 
method [25]. We have carried out the main part of the numerical analysis for the case of unit vortex number 
n = 1 in (8) as in Refs. [20] and [21]. In addition we have also studied more briefly, the case of n = 2. 

Starting with the case of vorticity n = 1, we have constructed monopole-antimonopole solutions for a large 
range of values of the coupling constant k. For vanishing coupling constant k the monopole-antimonopole solution 
corresponds to a non-Bogomol'nyi solution of the BPS system, for which our results are in good agreement with 
those of [21]. Our numerical analysis was carried out for the skyrmed model in the absence of the Higgs potential, 
namely with A = in (12). We did however check that the presence of nonvanishing A does not change the 
qualitative properties of our solutions. As expected the only effect it has is in the large r asymptotic region, 
where the modulus of the Higgs field for example, reaches its asymptotic value faster, namely exponentially. 

In Figure 1 we show the normalised energy of the solitons E/A-kt] and the energy Ei n f/4irr], of the monopole- 
antimonopole pair with infinite separation corresponding to twice the energy of a charge-1 monopole, as functions 
of the coupling constant k. As can be seen from Figure 1 the energy of the monopole-antimonopole solution is 
less than the energy of a monopole-antimonopole pair with infinite separation for all values of k. 




2 - 
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Figure 1: The energy of the monopole-antimonopole solution (solid line) and the energy of a 
monopole-antimonopole pair with infinite separation (dashed line), for n = 1. 

In Figure 2 we exhibit the modulus of the Higgs field \&(p, z)\ as a function of the coordinates p = \/ x 2 + y 2 and 
z for k = and n — 100. The zeros of |$(p, z)\ are located on the positive and negative z-axis at ±z ~ 2.1 for 
k = and at ±z w 1.5 for k — 100. The distance d of the two zeros of the Higgs field decreases monotonically 
with increasing k. 




Figure 2: The modulus of the Higgs field as a function of p and z for n = (a) and n — 100 (b), for n = 1 
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Asymptotically |$(p, z)\ approaches the value 1. But at the origin the value of the modulus of the Higgs field 
decreases monotonically with increasing n (see Figure 3). In the limit k — > oo \(j) \ w 0.015, and we expect the 
modulus of the Higgs field to be very small for \z\ < 4. 




2000 4000 6000 8000 10000 12000 14000 16000 18000 



Figure 3: The modulus of Higgs fields at the origin as a function of n, for n = 1 

In Fig ure 4 we show the energy density of the monopole-antimonopole solution as a function of the coordinates 
p = yjx 2 + y 2 and z for k — and k = 100. At the locations of the Higgs field the energy density posesses 
maxima. 

For small values of coupling constant k the equal energy density surfaces near the locations of the zeros of 
the Higgs field assume a shape close to a sphere, centered at the location of the respective zero (see Figure 4 
(a)). This presents further support for the conlcusion, that at the two zeros of the Higgs field a monopole and an 
antimonopole are located, which can be clearly distinguished from each other, and which together form a bound 
state. 



(a) 



(b) 





Figure 4: The dimensionless energy density as a function of p and z for k = (a) and k = 100 (b), for n = 1 

With increasing k the distance d between the monopole anti-monopole centres becomes smaller tending to a limit 
as k — ► oo. At the same time the spherical equal energy surfaces in Figure 4(b) become larger, and the equal 
energy density surfaces assume a shape that looks like the intersection of two spheres (see Figure 4 (b)), thus 
making it more difficult to distinguish the monopole from the anti-monopole. The dependence of the separation 
length d is given in Table 1 below as a function of k. 

Having exibited the qualitative properties of our dipole solutions, we give the values of the dipole moment 
that we calculated as a function of the coupling constant n, again in the Table 1. As expected, with decreasing 
d the dipole moment C m also decreases. 
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d 


4.19 


3.64 


3.49 


3.38 


3.29 


3.21 


3.16 


3.06 


2.54 


2.53 


Cm 


2.36 


2.27 


2.23 


2.19 


2.15 


2.11 


2.07 


2.02 


1.66 


1.65 



Table 1: Monopole anti-monopole separation d and dipole moment C m as a functions of k 

Finally we constructed solutions for the case of vorticity n = 2. Most of the qualitative properties of these 
solutions do not differ from those of the n — 1 case just described. The most noticable quantitative difference 
concerns the value of the modulus of the Higgs field at the origin, analogous with Figures 1(a) and 1(b). We do 
not exhibit here these analogous figures, but simply note that the the moduli of the Higgs fields at the origin are 
smaller than those in Figures 1(a) and 1(b) for the same values of the coupling constant k. 

Another difference, qualitative though expected, is that the surfaces of equal energy are not spheres centred 
on the z-axis but describe rings or tori around it. This is exhibited in Figures 5(a) and 5(b), analogously with 
Figures 4(a) and 4(b). 

(a) (b) 




Figure 5: The dimensionless energy density as a function of p and z for k = (a) and k = 100 (b), for n = 2 
Again, as k grows, the distiction between the monopole and anti-monopole rings gets blurred. 

5 Summary 

We have contructed axially symmetric solutions to a simple SU(2) skyrmed YM-Higgs model, with such boundary 
conditions that result in the description of a monopole anti-monopole pair with zero magnetic charge. These 
solutions have lower mass than two infinitely separated charge- 1 monopoles, and since they are characterised by 
zero magnetic charge, are not topologically stable. 

When the usual boundary conditions are imposed, the skyrmed SU{2) YM-Higgs model employed here sup- 
ports mutually attractive monopoles, including axially symmetric charge-2 monopoles. This is in contrast to the 
Georgi-Glashow model studied in [21] where due to the Higgs potential the monopoles are mutually repulsive [12]. 
Nevertheless, the qualitative features of the monopole anti-monopole solutions in the two models arc similar. In- 
creasing the Skyrmc coupling constant k in the present model results in the approaching of the monopole and 
the anti-monopole centres down to a limiting value 2.53 as k — > oo, just as it does to the limiting value 3.0 as 
A — > oo in the Georgi-Glashow model A being the Higgs coupling constant. (Our results are for A = 0.) 

Another parallel property in the two models is the changing dipole moment with respect to the change in 
the Skyrme coupling constant k and the Higgs coupling constant A, in the two models respectively. Specifically 
in the present model the magnetic moment decreases with increasing k, with limiting value 1.64, while in the 
Georgi-Glashow model it decreases with increasing A, with limiting value 1.55, in the same units. 

Finally, we studied also the case of a zero charge monopole which has vortex number n = 2 rather than n = 1. 
The qualitative properties again stay unchanged. The most noticable quantitative difference of the n = 2 soltion 
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is that the value of the modulus of the Higgs field at the origin is smaller than that of the n = 1 solution, for the 
same value of k, and, the distance between the two centres is also smaller. For example at k — 25 the distance 
d = 3.38 for the n = 1 solutions while that for the n — 2 is d = 1.33. 
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